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Spectral signatures of thermal spin disorder and excess Mn in half-metallic NiMnSb 
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Effects of thermal spin disorder and excess Mn on the electronic spectrum of half-metallic NiMnSb 
are studied using first-principles calculations. Temperature-dependent spin disorder, introduced 
within the vector disordered local moment model, causes the valence band at the T point to broaden 
and shift upwards, crossing the Fermi level and thereby closing the half-metallic gap above room 
temperature. The spectroscopic signatures of excess Mn on the Ni, Sb, and empty sites (Muni, 
Mnsb, and Mue) are analyzed. Muni is spectroscopically invisible. The relatively weak coupling of 
Mnsb and Mue spins to the host strongly deviates from the Heisenberg model, and the spin of Mue 
is canted in the ground state. While the half-metallic gap is preserved in the collinear ground state 
of Mnsb, thermal spin disorder of the weakly coupled Mnsb spins destroys it at low temperatures. 
This property of Mnsb rnay be the source of the observed low-temperature transport anomalies. 


Half-metallic ferromagnets urn with a band gap in 
one spin channel can serve as sources of highly spin- 
polarized currents in spintronic devices [5l [6]. While 
many materials, including a number of Heusler com¬ 
pounds mm, are predicted by band structure calcula¬ 
tions to be half-metallic, experimental verification is com¬ 
plicated by defects, thermal spin excitations, and surface 
non-stoichiometry [7]. Correlation effects can also lead 
to spin depolarization already at zero temperature [8]. 

If a perfect crystal is half-metallic at T = 0, there may 
be a region of low T and weak disorder where transport 
is effectively single-channel [6]. However, a crossover to 
two-channel [9] conduction should occur if a quasiparti¬ 
cle band crosses the Fermi level Ep Ai the gapped spin 
channel. Thermally-induced incoherent spectral density 
in the gap can also lead to transport anomalies. 

A bulk-sensitive positron annihilation measurement 
m found NiMnSb to be half-metallic to experimental 
accuracy, but other measurements paint a more compli¬ 
cated picture. An anomaly in the transport properties 
was observed near 80 K mm- It was suggested [13] 
that it could be due to the longitudinal fluctuations of 
the Ni spin moments, but the effect was overestimated 
by improper phase space integration El- Recent stud¬ 
ies indicate that rather than being intrinsic, the 80 K 
anomaly is only found in samples with excess manganese 
na nnj. This suggests the anomaly is not related to 
the non-quasiparticle states [8], which also lack a clear 
parameter of order 10 meV associated with them. It 
was suggested that thermal excitations on the interstitial 
Mue atoms may be involved [16] • Here we study the in¬ 
fluence of thermal spin fluctuations and excess Mn on the 
single-particle electronic spectrum of NiMnSb and iden¬ 
tify Mnsb as the most likely source of the 80 K anomaly. 

For averaging over spin disorder in electronic struc¬ 
ture calculations, we used the adiabatic disordered lo¬ 
cal moment (DLM) method [18], which provides excel¬ 
lent agreement with supercell averaging for the density 
of states (DOS) and electric resistivity of ferromagnetic 


metals [iniisn]. The effects of spin disorder at low tem¬ 
peratures can not be described by the binary pseudo-alloy 
approximation with collinear “up” and “down” atomic 
spin moments [21], which is justified only for the param¬ 
agnetic state [18]. Instead, we employ proper averaging 
[22] over all directions of the Mn spins. Below we refer 
to this approach as the vector DLM (VDLM) model. 

In ideal NiMnSb we used the single-site mean-field dis¬ 
tribution function p{0) oc ex.p{acos0) for Mn spin direc¬ 
tions. The parameter a = where Rw is the 

Weiss field and p the local spin moment, is mapped to the 
temperature using the experimental M{T) curve [1711^ . 
Both spin disorder and excess Mn defects are treated 
within the coherent potential approximation (CPA) ap¬ 
plied within the Green’s function-based linear muffin-tin 
orbital (GF-LMTO) method [24H27] . The azimuthal an¬ 
gle (f in the CPA equations is integrated out analytically, 
while the 0 dependence is discretized using the Gauss- 
Legendre quadrature. The potentials for all atoms (and 
0 angles) are determined self-consistently by embedding 
the GPA self-consistency loop into the density-functional 
iteration. To enforce magnetic self-consistency, a trans¬ 
verse constraining field [28] is introduced for each 0. The 
experimental lattice constant of 5.92 A and the general¬ 
ized gradient approximation m are used throughout. 

The electronic structure in the presence of disorder is 
represented by the energy- and k-resolved GPA spec¬ 
tral function [25]. It (and DOS) is resolved by spin 
in the global reference frame and is accessible through 
spectroscopic measurements like spin-polarized angular- 
resolved photoemission (ARPES). Our focus will be on 
the minority-spin spectral function, which is gapped at 
T = 0 in the ideal crystal. A small 10“^ Ry imaginary 
part is added to energy to improve GPA convergence; its 
effect is negligible compared with disorder broadening. 

First we consider the spectroscopic effects of thermal 
spin disorder in defect-free NiMnSb. Fig.j^a) shows the 
minority-spin band structure at T = 0, and Fig.[^b) the 
spectral function at an elevated temperature correspond- 
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ing to the reduced magnetization M(T)/M(0) = 0.80. 

Mixing with majority-spin bands is apparent at 1.2-1.6 
eV below where those bands and avoided crossings 
show up in Fig. [^b). The bands that tend to acquire 
a significant admixture of the opposite spin character at 
T 7^ 0 are the flat bands with small dispersion. We do 
not, however, see any trace of the dispersive majority- 
spin bands that cross the Fermi level. The spectral weight 
is, of course, transferred from those bands to the minor¬ 
ity spin, but this transferred weight is incoherent and 
spreads over the entire bandwidth. Bloch electrons in 
dispersive bands behave according to the Stoner picture 
and “observe” an averaged effective field, because adjust¬ 
ing their wavefunctions locally to the directions of indi¬ 
vidual spins costs a lot of kinetic energy. This general 
property is favorable for spintronic applications, because 
the quasiparticles near Ep retain their spin polarization 
in the presence of spin disorder. 

Fig. [TJb) shows that the valence-band maximum 
(VBM) triplet at F is strongly broadened and shifted 
upward by spin disorder. The conduction band mini¬ 
mum (CBM) at X is also broadened and shifted upward 
but at a smaller rate. Fig. plots the shift and broad¬ 
ening of VBM and CBM as a function of T. The VBM 
comes within its half-width of Ep at about 400 K. This is 
where strong transport spin depolarization and crossover 
to conventional two-channel conduction can be expected. 
Although this temperature depends on the initial posi¬ 
tion of at T = 0, it appears that spin disorder can 
not lead to low-temperature anomalies in ideal NiMnSb. 

The spin disorder-induced upward shift of the 
minority-spin VBM states is contrary to the expecta¬ 
tion based on the simple exchange splitting picture. This 
counterintuitive band shift can be understood by exam¬ 
ining the hybridization that forms this band in Fig.j^c). 
As seen in Fig.j^a) (and confirmed by full-potential cal¬ 
culations), the VBM triplet at F is predominantly com¬ 
posed of Mn t 2 g states with an admixture of Sb bp and 
Ni t 2 g states. If the 3d states of Mn and Ni are removed 
from the basis set, the remaining bands deriving primar¬ 
ily from Sb bp states (green in Fig. [^c)) extend from —6 
to 2 eV. The F point is antibonding for these bands. If, 
on the other hand, all Sb states except bs are removed 
from the basis set, we end up with the bands for which 
the Mn character is shown by red {t 2 g) and blue (e^) 
color in Fig.j^c). Comparing with Fig.j^a), we see that 
the hybridization of the Mn t 2 g and Sb 5p-derived states 
[1] opens the gap and forms the VBM at F. Spin disor¬ 
der reduces this hybridization by adding the spin overlap 
factors to the hopping matrix elements, resulting in the 
upward shift of the VBM triplet seen in Fig. [^b). 

As the 80 K anomaly may be due to excess Mn usmsi, 
we now consider the effects of excess Mn on the Ni, Sb, 
and empty sites (Mn^i, Mnsb, MnE). Various defects in 
NiMnSb were studied using supercell calculations |30l[3l] 
and CPA |3T1 122] . While the energies of atom-pair swaps 


m are related to the formation energies of different de¬ 
fects, the concentration of each non-stoichiometric defect 
depends on the chemical potentials of the relevant ele¬ 
ments and thus on the material synthesis protocol. A 
polarized neutron diffraction measurement [33] of a low- 
quality crystal found a large concentration of Mnsb and 
SbMn defects. Since the Mn-Sb swap energy is rather 
high m, the formation of these defects could be due to 
nonequilibrium growth conditions. 

We assume excess Mn is randomly distributed over the 
given sublattice and treat each type of defect separately 
within CPA. First we used VASP [34| to relax the atomic 
structure for a 48-atom (2x2x1) host supercell [35] with 
one excess Mn defect and examine its magnetic coupling 
to the bulk. For Mn^i the antiparallel spin alignment 
is very stable in agreement with Ref. |36| while the self- 
consistent solution for the parallel alignment could not be 
obtained. Fig.j^shows the total energy E{0) as a function 
of the angle 0 made by the spin of MnE or Mnsb with the 
magnetization. In both cases there are strong deviations 
from the Heisenberg model predicting E{0) oc cos^. For 
MnE we confirm the earlier result m that E{7r) is lower 
than E{0). However, the energy minimum is reached for 
a canted spin at 0 ^ 130°. There is a high-to-low spin 
crossover near 0 = it/A. For Mnsb the minimum is at 
^ = 0, but the range of energy variation is only 70 meV. 

Since the magnetic coupling of MnE and Mnsb spins to 
the bulk is weak, their disordering can lead to anomalies 
at relatively low temperatures. We now use CPA to study 
the spectral properties of these defects. We start with the 
6.25% defect concentration as in the VASP supercells, 
which is also typical for low-quality samples [33] . 

First we performed an LMTO calculation for each 
supercell, taking the optimized atomic positions from 
VASP and setting the sphere radius of excess Mn to the 
same value as the substituted atom The partial 

DOS (PDOS) in these calculations are similar to VASP 
results, and they give us the correct charge and spin mo¬ 
ments in the atomic sphere for the excess Mn atom. In 
CPA calculations with the ideal lattice, a screening cor¬ 
rection to the Madelung potential is applied with a 
coefficient that is adjusted to provide the correct charge 
on Mn. In addition, the local part of the effective GGA 
magnetic field for the MnE defect was scaled by a factor 
1.2 to obtain good agreement with the correct local spin 
moment for MnE at both ^ = 0 and tt. This correction 
compensates for the lack of atomic relaxations in CPA. 
The PDOS with these corrections agree well with super¬ 
cell calculations up to disorder broadening (see Fig. [^for 
Mnsb at ^ = 0; the agreement is also good for 0 = tt). 

The spectral function with Mn^i defects (not shown) 
is only slightly modified compared to ideal NiMnSb (Fig. 
0a)), and the broadening of the bands is very small. 
Since the spin of Mn^i is strongly coupled to the host 
magnetization [36], this defect is spectroscopically invis¬ 
ible and can not lead to low-temperature anomalies. 
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FIG. 1. (a-b) Minority-spin spectral functions in NiMnSb. (a) T = 0. (b) Spin disorder (VDLM) at M(T)/M(0) = 0.80. 
The coloring reflects the spectral weight contributions of Mn (red), Ni (blue), and Sb (green, barely seen), (c) Hybridization 
analysis: Red and blue coloring show, respectively, Mn t 2 g and Mn Cg orbital character when all Sb states except 5s are removed 
from the basis set. Green coloring shows the Sb character when the 3d states of Mn and Ni are removed from the basis set. 
(All charges are taken from the calculation with the full basis set.) Energy in all panels is in eV measured from Ep- 



FIG. 2. Temperature dependence of the band center and half¬ 
width of the minority-spin bands at T and X points at the 
edges of the half-metallic gap in NiMnSb. 



FIG. 3. Energy of a supercell with one Mue or Mnsb defect 
as a function of its spin orientation 0 relative to the host 
magnetization (from VASP). HS, LS: high spin, low spin. 


Spectral functions for Mue at ^ = tt and in the fully 
spin-disordered (SD) state are shown in Fig.a-b). For 
the SD state the 0 integration in CPA is combined with 
full charge self-consistency. The spin moment of MnE 
varies from 0.84 to 3.04 fiB between the points closest 
to 0 and tt in good agreement with VASP calculations. 
PDOS at 6> = 0, TT, and SD state are shown in Fig. 

Fig. a-b) and show that spin disorder on MnE in¬ 
troduces a substantial minority-spin spectral weight 
inside the half-metallic gap. Its magnitude at E'e is about 
13% of the majority-spin DOS N^. However, if MnE is 



E-Ep, eV 

FIG. 4. Partial DOS for Mnsb in NiMnSb with a 6.25% 
occupation of the Sb site by Mn, parallel spin alignment. 
Solid (blue) lines: GPA; dashed (red) lines: VASP. Black dot- 
dashed line: total minority-spin DOS per formula unit (iV^), 
disordered spin of Mnsb- Inset: at 40 K, 60 K, 80 K, 120 

K, and 160 K as labeled (same units on axes). 


already canted at T = 0 as our calculations suggest (Fig. 
1^, it is likely not the source of the low-T anomaly. 

Fig. i c-d) shows the spectral functions for Mnsb at 
^ = 0 (ground state) and in the SD state. The spin 
moment of Mnsb varies from 3.82 to 4.18 {Ib between 
angles closest to 0 and tt. Fig. [^c) shows, in agreement 
with Ref. [30l that Mnsb preserves the half-metallic gap in 
the ground state, although the VBM moves closer to Ep 
compared to ideal NiMnSb. Further, in Fig. id) and 
1^ (dot-dashed line and inset) we see that spin disorder 
introduces a large spectral density inside the gap near F 
above the VBM, which is much larger compared to Mue 
(F ig. [^b)). These features are due to the replacement of 
Mn-Sb bonds by Mn-Mn; recall that Mn-Sb hybridization 
pushes the VBM states downward at F (Fig. Ba,c)). 

The inset in Fig. shows near Ep for several tem¬ 
peratures, which were obtained using the VDLM model 
with the spin distribution function p{0) oc exp[—E(^)/T] 
for Mnsb with E{0) from VASP (Fig. [^. We see that 
thermal spin disorder on Mnsb generates a consider- 
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FIG. 5. Minority-spin spectral function in NiMnSb with 6.25% of: (a-b) Mue, (c-d) Mnsb- The spins of Mue are: (a) 
antiparallel to magnetization; (b) randomly oriented. The spins of Mnsb are: (c) parallel to magnetization; (d) randomly 
oriented. Coloring is similar to Fig. [^a); Mue added to Mn (red), Mnsb to Sb (green). Energy is in eV measured from Ep- 



FIG. 6. Total DOS for NiMnSb with a 6.25% occupation of 
the empty site by Mn (CPA). FM (AFM): Mue is parallel 
(antiparallel) to the magnetization. VDLM: fully disordered 
spin on Mue. Upper inset: magnified region near Ep- Lower 
inset: partial DOS for Mue. (Same units on axes in insets.) 

able at Ep at rather low temperatures; for example, 
^ 0.26A^^ at 80 K. The spectral spin polarization at 
Ep, P = N^)/ + A^ 4 ,), is mapped in Fig. as 

a function of temperature and Mnsb concentration. It 
confirms that the quick drop of P occurs in the 70-100 
K range near 5% Mnsb, consistent with observations of 
Ref. [151 At lower concentrations the crossover is shifted 
to higher temperatures; this is a testable feature of the 


thermal depolarization mechanism associated with Mnsb- 



FIG. 7. Spectral spin polarization P at Uf as a function of 
Mnsb concentration x and temperature. 

In conclusion, we have studied the influence of thermal 
spin disorder and excess Mn on the electronic structure of 
NiMnSb. Thermal spin disorder broadens and shifts the 
minority-spin VBM at F upwards due to the unmixing 
of Mn t 2 g states from Sb states. While Mn^i is spec¬ 
troscopically invisible, Mue and Mnsb have strong spec¬ 
tral signatures. Mnsb is identified as a possible source of 
the low-T anomaly: its spin couples weakly to the bulk, 
and easily excited spin disorder strongly contaminates 
the half-metallic gap near F. These predictions can be 
tested through a combination of transport and spectro¬ 
scopic measurements. 
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